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\Q . We prove the finite-time collapse of a system of N classical fields, which are described by N coupled 

nonlinear Schrodinger equations. We derive the conditions under which all of the fields experiences 
this finite-time collapse. Finally, for two-dimensional systems, we derive constraints on the number 

CI of particles associated with each field that are necessary to prevent collapse. 

O ' 

' PACS numbers: 

Q: 

The nonlinear Schrodinger equation (NLSE) arises in many diverse areas of physics, such as nonlinear optics, gravity 
waves in deep water and, recently, trapped dilute Bose-Einstein gases where it has been important in explaining many 
' ?H '" °f the static and dynamic properties of dilute Bose-Einstein condensates (BECs) [H, 0- It has been proven that, 
under certain conditions, the NLSE exhibits finite-time collapse, which is the occurrence of a finite time singularity 
^3 ' in the field. In this paper, we generalize this proof to show that finite-time collapse occurs in N-coupled nonlinear 
O | Schrodinger equations and we identify the parameter regime and initial conditions required for this collapse. We 
motivate the work of this paper by referring to the state-of-the-art in experimental research on BECs. 

The NLSE was shown, in the context of Langmuir waves in plasma physics, to undergo finite-time collapse for 
systems with a negative coupling constant Q. Finite-time collapse was indicated by the absence of variance in the 
| field described by the NLSE at a finite time. Proof of this finite-time collapse of the "focusing NLSE" has been 
Ch ' important for understanding the stability of trapped single-component BECs. In particular, this proof implies that a 
condensate formed by atoms with a negative coupling constant, i.e. atoms that attract each other, would in theory 
develop a singular density; in reality, three-body collisions disallow the development of an infinite density (In 
experimental setups involving a small number of attractive atoms, however, an external potential can be used to 

■ stabilize the dilute BEC system even though the coupling constant is negative.) 

We move now to the generalization of this proof of the focusing NLSE's finite-time collapse to that of N coupled 
NLSEs describing a system of coupled classical fields. In generalizing the proof, we shall also determine the range 
\Q ' of initial conditions and the parameter regime over which all of these coupled fields exhibit instability. The present 
exercise is of particular importance given the shift in focus of BEC experiments from single-component condensates 
to condensates of miscible mixtures. Current experiments are able to achieve two component condensates (see for 

■ example 5, 6]), and the prospects for realizing condensates of many more components are good given, for instance, 
the recent achievement of condensed Ytterbium Q, which has five bosonic isotopes. 

Coupled NLSEs have been successfully employed to model many dynamic and static properties of binary mixtures 
in recent dilute Bose-Einstein condensate experiments (see for example @, Q). As a generalization of this approach 
I ■ for multicomponent systems, let us consider the following N coupled NLSEs describing the behavior of the classical 
'"O ' fields tpi of the system: 
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where gij are the coupling constants (which, in the BEC context, can be written in terms of the relevant two-body 
jj^ ■ scattering length). For simplicity we will assume g^ — gji and no external potentials. 

To find the conditions that result in finite-time collapse of all classical fields, we examine the sum of the variances 
of the fields in the system, which can be expressed as 

N N . 

V(t) = X>i(t) = X; / drr^(r,t)\ 2 . (2) 

1=1 i—l 

When V(t) — 0, since the variance of each given field, Vi(t) is necessarily positive, all of the coupled classical fields 
must have a variance of zero, i.e. all fields become singular. 

To prove that V(t) — > in a finite time, using eq. [TJ we calculate the variance for a given field (through successive 
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integration by parts) 



d^V^t) = / dr 
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where, for simplicity, we have ignored the angular variables and, for completeness, we have kept the dimension of 
space d as a variable. Simplifying further using integration by parts, we arrive at the crucial result (which was not at 
all obvious for the multicomponent system) for the proof of the finite time collapse for coupled NLSEs: 

d 2 V(t) = 8H + 4(d- 2)U(t). (4) 

Here the Hamiltonian, which is a constant of motion of this system, is given by 

H = T(t) + U{t) (5) 

where T(t) corresponds to the kinetic energy of the system and is given by 

N 

T «=E M 3 r|V^M)| 2 ; (6) 

i=l J 

and U (t) corresponds to the potential energy of the system and is given by 

U(t) = \J d 3 r^(r,t)| 2 |^(r,t)| 2 , (7) 

summing over repeated indices. 

Assuming d = 3, we must impose the condition U(t = 0) < in order to prove finite-time collapse using the variance 
method presented in this paper. This implies that the matrix g^j must be negative semi-definite (llj . If d = 2, there 
is no need to assume U(t) is negative. Then, following the logic presented in [9| for the single NLSE, we impose on 
eq. ((4]) the condition that gij be negative semi-definite for d — 3 (no such condition is imposed for d — 2) to arrive at 

2 V(t) < 8H. (8) 

Since H is a constant of motion, we can integrate this equation to obtain 

V(t) < AHt 2 + V'(0)t + V(0). (9) 

The right-hand side will go to zero in a finite time - which, as argued above, implies that every one of the fields ipi 
must become singular within a finite time - if any one of the following three sets of conditions is true: 

1. H < 

2. H = and V'(0) < 



3. H > and V'(0) < -iy/HV(0) 

Note that the time given in eq. © is only an upper bound on the collapse time of the entire system. 

The fact that the final term in eq. ((3]) is proportional to U(t), which was not expected a priori, greatly simplifies 
the proof of finite time collapse of classical fields whose behavior is determined by coupled NLSEs. This leads to 
the simple conclusion that if the potential energy, U(t), is sufficiently negative to overcome the kinetic energy, T(t), 
collapse will occur. 

As an example, let us assume initially that every one of the fields is uniform, i.e. T{t = 0) = 0. In this case, a 
sufficient condition for finite-time collapse of all fields is that the matrix be negative semi-definite, which assures 
U(0) < 0. Therefore, as the sum of T(0) and U(0), H(0) < 0, which satisfies the first set of conditions above. 

If d = 2, one can derive strict contraints on the number of particles in each component given by iVj(i) = / dr\ijji(r, t)\ 2 
in order for the coupled classical fields not to collapse. For clarity, let us first derive the analogous constraint on the 
number of particles for the collapse of a single field obeying the NLSE. 

The proof relies on the Gagliardo-Nirenberg inequality given by (note that the discussion below assumes d = 2) 

di#M)| 4 < CN(t)T(t) (10) 
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where C is a positive constant determined by the appropriate Euler-Lagrange equations [10(. Note that we have 
temporarily removed the subscript as we are dealing here with only one non-coupled component. Using this inequality, 
one can show that 

H > J dr\iP(r,t)\ 4 

Therefore, if g > 0, or if N(t) < —2/gC and g < 0, then the potential energy U{t) is bounded by a finite amount, 
which implies that the kinetic energy is also bounded, and hence collapse is not possible. 

One can generalize this argument to the case considered in this paper of N classical fields governed by coupled 
NLSEs given by eq. ((T|) for d = 2. Using the inequality given by eq. (fltj|). one can show that 

H(t) > J drA y #,(r,i)| 2 |^(r,;)| 2 (12) 

where A y - = <5y /CNi(t) + gij/2 and <5 y - = 1 when i = j and zero otherwise. If Ay is positive definite [l]|, then there 
will be no collapse as the potential energy is bounded by a finite amount. This condition, assuming <7 y is given, puts 
constraints on the number of particles in each coupled field to ensure there is no collapse in d = 2, analogous to the 
single-component case considered above. Note that if Ay is not positive definite, then this analysis does not give 
information about whether or not collapse will occur. 
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